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ABSTRACT 
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This  report  outlines  the  following  four  areas  of  interest: 

1.  The  concept  of  the  semiadiabatic  engine, 

2.  The  creation  of  a  mathematical  model  to  establish  heat  transfer  and 
temperature  profiles  for  such  engines, 

3.  The  establishment  of  operational  boundary  conditions  for  a  semiadiabatic 
engine,  and 

*T.  Discuss  the  results  of  an  endurance  test  of  ceramic  coated  engine  parts. 
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INTRODUCTION 


V 

The  following  document  is  a  final  report  for  Grant  Number  DAAG29- 
79-G-0021  funded  by  the  U.S.  Army  Research  Office.  This  project  was 
started  1  January  1979  and  had  a  duration  of  two  years. 

In  general,  the  project  was  1)  to  develop  a  mathematical  model 
describing  heat  transfer  and  temperature  characteristics  in  ceramic 
coated  engine  components,  2)  measure  operational  temperatures  for 
such  parts  with  various  coating  thickness,  and  3)  continue  a  previous 
endurance  test  of  ceramic  parts. 

/V 

THEORETICAL  BACKGROUND 


PURPOSE 

Since  the  time  of  Rudolph  Diesel,  engineers  have  been  trying  to 
develop  engines  that  waste  less  energy, and  thereby  develop  more  power 
with  a  given  quantity  of  fuel.  It  is  widely  accepted  that  heat  trans¬ 
ferred  to  the  cooling  system  of  an  engine  represents  wasted  energy. 

On  the  other  hand,  the  only  practical  way  known  to  control  engine  temp¬ 
eratures  is  to  provide  such  cooling.  In  reciprocating  engines,  the 
design  of  pistons  and  valves  has  always  required  careful  consideration 
of  temperature  as  related  to  component  strength  and  to  lubrication. 

But  the  cooling  necessary  to  keep  these  components  from  overheating  not 
only  reduces  engine  power  production,  but  consumes  power  to  drive  the 
cooling  mechanism. 

It  was  recognized  by  the  author  that  an  insulating  surface  on  com¬ 
bustion  chamber  surfaces  offered  an  approach  to  the  cooling  problem  by 
reducing  heat  transfer  to  the  cooling  system  and,  at  the  same  time, 
provided  lower  temperatures  of  metallic  components.  Lower  metallic 
temperatures  allow  redesigning  of  components  for  less  weight  and  lower 
frictional  drag. 

From  a  design  standpoint,  therefore,  lower  piston  temperatures  en¬ 
able  lower  piston  weight,  which  leads  to  reduced  inertia  loadings  thus 
allowing  smaller  and  lighter  components  such  as  crankshaft  bearings. 
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connecting  rods,  and  crankcase  structure.  The  lower  temperatures  also 
allow  a  new  freedom  in  selecting  lubricating  oils  that  can  reduce  vis¬ 
cous  drag  and  windage  losses  within  the  crankcase. 

A  second  significant  benefit  that  may  result  from  ceramic  coatings 
is  the  reduction  of  the  failure  rate  in  valves  and  pistons,  and  there¬ 
fore  an  increase  in  engine  reliability.  While  replacement  costs  for 
individual  engine  parts  are  not  unreasonable,  the  time  of  failure  and 
the  labor  required  are  often  critical  factors  in  vulnerable  situations. 

Engine  failure  is  particularly  critical  in  military  applications. 
Often  military  engines  are  destroyed  by  a  combination  of  operational 
factors  in  or  near  a  combat  zone.  For  example,  it  is  not  uncommon  for 
the  operator  to  refuel  with  any  available  fuel  and  then  operate  the 
vehicle  at  peak  power  for  extended  time  intervals.  Either  a  spark  ig¬ 
nition  (S-I)  or  a  compression  ignition  (C-I)  engine  operated  in  this  mode 
will  develop  severe  knock.  Knock,  in  turn,  dramatically  increases  con¬ 
vection  boundary  layer  heat  transfer  which  may  raise  valve  and  piston 
temperature  to  a  point  of  failure. 

Figure  1,  following,  is  a  photograph  illustrating  a  diesel  exhaust 
valve  failure  caused  by  excessive  thermal  stress  (hoop  stress  failure). 

A  more  common  failure  occurs  when  the  valve  head  separates  from  the  stem 
("valve  swallowing"),  after  high  temperature  stress  and  erosion.  This 
second  failure  mode  usually  occurs  at  high  power  and  high  rpm  conditions 
resulting  in  substantial,  if  not  complete,  engine  destruction. 

Two  examples  of  piston  damage  resulting  from  high  temperature  operation 
are  shown  in  figures  2  and  3.  Both  are  caused  by  knock  induced  heat 
transfer.  The  second,  however,  was  aggravated  by  preignition  at  the 
spark  plug  center  electrode. 

Insulated  coatings  are  expected  not  only  to  economically  eliminate  all 
of  the  above  types  of  engine  failure,  but  also,  to  increase  engine 
efficiency  and  reduce  overall  engine  weight. 

Based  on  the  information  above,  it  appears  that  one  small  breakthrough 
in  controlling  heat  transfer  offers  the  opportunity  for  major  engine 
improvement.  For  these  reasons,  it  is  highly  desirable  to  gain  a  better 
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understanding  of  how  heat  is  transferred  in  ceramic  coated  internal 
combustion  engine  parts. 

PREVIOUS  WORK 

Shortly  after  the  principal  investigator  realized  the  possible 
benefits  of  insulated  combustion  chamber  surfaces,  he  wrote  a  math¬ 
ematical  model  to  investigate  coated  engine  performance.  The  model 
did  predict  both  increased  efficiency  and  lower  component  temperatures. 
Following  analysis  of  the  model,  he  participated  in  two  engine  test 
programs  involving  coated  parts. 

The  first  engine  test  program  was  a  multicylinder  test  at  the 
Johnson  Space  Center  in  Houston.  The  major  outcome  of  this  test  was 
proof  that  plasma  sprayed  ceramic  coatings  could  endure  the  torturous 
environment  of  an  engine  combustion  chamber  for  extensive  periods  of 
time  without  significant  deterioration. 

The  second  test  program  was  a  U.S.  Army  funded  study  at  Oklahoma 
State  University  which  investigated  engine  performance  characteristics. 
Both  spark  ignition  (S-I)  and  compression  ignition  (C-I)  performance 
were  investigated  at  various  compression  and  equivalence  ratios.  These 
tests  incorporate  a  thin  coating  (0.025")  of  yttrium  stabilized  zircon¬ 
ium  oxide  on  the  piston,  valves,  and  cylinder  head  surfaces  of  a  Cooper¬ 
ative  Fuel  Research  (CFR)  engine.  In  addition,  tests  were  undertaken 
to  investigate  how  various  catalysts  applied  to  these  coatings  will 
affect  combustion.  The  details  of  this  project  are  found  in  the  final 
report, project  No.  P-14099-E  entitled  "Improved  Engine  Performance  Through 
Heat  Transfer  Control",  available  by  contacting  the  U.S.  Army  Research 
Office,  Research  Triangle,  North  Carolina. 

CERAMIC  COATINGS 

Several  ceramic  substances  have  been  developed  in  the  last  few  years 
that  exhibit  excellent  insulating  properties.  These  substances  can  be 
sprayed  on  metal  surfaces  by  flame  or  plasma  methods  and  also  can  be 
applied  by  a  device  known  as  a  detonation  gun.  Some  properties  of  various 
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coatings  and  metals  follow  in  table  1. 

To  reduce  heat  transfer  in  an  engine  cylinder,  two  physical 
coating  properties  are  of  prime  importance.  First,  the  material 
must  adhere  to  the  metallic  subsurface  without  any  significant  in¬ 
terface  stress  over  a  wide  range  of  temperature.  This  necessitates 
that  the  coefficient  of  thermal  expansion  of  the  coating  should  be 
of  the  same  order  of  magnitude  as  the  base  structure.  Second,  it 
should  have  as  low  a  thermal  conductivity  as  possible  in  order  to 
effectively  control  heat  transfer.  Properties  of  lesser  importance 
are  specific  heat,  compressive  strength,  and  adherence.  Of  all  the 
coating  materials  investigated, yttrium  stabil ized  zirconium  oxide 
best  fits  the  requirements  for  engine  combustion  chamber  insulation. 

Adherence  of  plasma  and  flame  sprayed  zirconium  oxide  to  cast 
iron  and  aluminum  surfaces  is,  however,  in  some  cases  not  good,  and 
coating  separations  will  occur  if  proper  techniques  are  not  observed. 
Metallic  surfaces  must  be  clean  and  free  of  oil  prior  to  coating.  In 
some  cases.it  may  be  necessary  to  preheat  cast  iron  surfaces  to  expel 
residual  oil  from  machining  process.  To  promote  good  bonding,  it  is 
suggested  that  metallic  surfaces  be  grit  blasted  and  given  a  flash  coat 
of  nickel  aluminate  before  application  of  the  zirconium  oxide  coating. 


PROJECT  SCOPE 

The  preceding  section  has  outlined  the  historical  use  of  ceramic 
coatings  to  improve  engine  performance  and  to  make  engines  more  toler¬ 
ant  of  low  grade  fuel.  This  project  was  undertaken  to  establish  a  data 
base  to  assist  design  engineers  in  the  adoption  of  ceramic  coatings. 

To  achieve  this  goal,  a  two  year  project  was  implemented  to: 

1)  create  a  mathematical  model  that  can  predict  piston 
and  exhaust  valve  heat  transfer  as  a  function  of 
engine  design  parameters  and  coating  characteristics, 

2)  establish  operational  temperature  profiles  for 
coated  pistons  and  valves,  and 


3)  an  additional  assignment  for  this  project  was  to 
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Table  1 

PHYSICAL  PROPERTIES 


Material 

Thermal 

Conductivity 

Coefficient 

of 

Expansion 

Bulk 

Density 

Specific 

Heat 

Melting 

Temperature 

BTU 

in  x  10"6 

gm 

BTU 

°F 

hr  ft  °F 

in  °FX  IU 

cm3 

#  °F 

i 

1 

1 

png 

wmm 

(°C) 

■RilSSiiBI? 

■  ■ 

■ 

^Zirconium 

Oxide 

0.67 

(2.77  x  10"3) 

5.4 

(3.0) 

5.2 

B 

4500 

(2482) 

Aluminum 
j  Oxide 

i 

1.58 

(0.01) 

4.1 

(2.28) 

3.3 

■ 

3600 

(1982) 

Chrome 
|  Oxide 

| 

5.0 

(2.78) 

n 

■ 

3000 

(1649) 

Aluminum 

1 _ _ _ 

■■ 

13.0 

(7.22) 

D 

. 

mi 

Cast  iron 

27 

(0.11) 

6.5 

(3.61) 

7.88 

0.11 

(0.011) 

2800 

(1538) 

Source :"Rokide  Ceramic  Spray 
Coating ", 

The  Norton  Company,  1961. 
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continue  an  endurance  test  for  ceramic  coated 
engine  parts.  Prior  to  the  start  of  this  pro¬ 
ject  an  engine  had  undergone  1000  hours  of 
operation  simulating  highway  conditions.  A 
second  1000  hours  of  operation  was  established 
as  a  third  goal  of  this  project. 


THE  MATHEMATICAL  MODEL 


GENERAL 

Engine  pistons  and  exhaust  valves  are  three  dimensional  objects 
having  configurations  that  can  typically  be  portrayed  by  figures  4 
and  5,  respectively. 

The  piston  and  valve  combustion  chamber  surfaces  receive  a 
cyclic  heat  input  during  each  combustion,  expansion  and  exhaust  process. 
This  heat  flow  depends  primarily  on  the  power  loading  of  the  engine 
and  the  combustion  characteristics.  These  same  surfaces  are  also  sub¬ 
jected  to  a  cyclic  cooling  during  the  intake  and  compression  processes. 
The  cooling  tendency  varies  mainly  with  air-fuel  ratio  and  engine  power 
level . 

Primary  piston  cooling  is  through  the  piston  rings  and  heat  trans¬ 
fer  to  lubricating  oil,  both  being  somewhat  cyclic  in  nature. 

Exhaust  valve  heat  rejection  occurs  primarily  to  the  valve  seat  and 
valve  guide.  Both  of  these  sinks  change  characteristics  with  time  and 
therefore  are  transient  in  nature. 

While  both  the  piston  and  exhaust  valve  are  transient  heat  transfer 
devices,  they  may  be  treated  as  steady  state  for  practical  purposes. 

The  insulating  effect  of  a  ceramic  surface  tends  to  isolate  the  core 
metal  structure  from  temperature  fluctuations,  in  other  words,  it  is 
believed  that  the  interior  temperature  of  both  such  pistons  and  valves 
is  controlled  more  by  the  constant  temperature  sink  and  the  physical 
properties  of  the  metal  than  by  the  cyclic  heat  load.  This  theory  was 
verified  in  the  laboratory;  details  and  test  data  confirming  this  will 
be  presented  later. 

As  seen  in  figure  5,  the  exhaust  valve  is  symmetrical  and  can, 
therefore,  be  treated  as  a  two  dimensional  heat  transfer  body. 
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The  piston  is  also  basically  symmetrical  with  the  exception  of 
the  wrist  pin  boss.  This  boss,  however,  is  typically  located  well 
below  the  ring  area  and,  therefore,  is  usually  considered  to  have 
little  contribution  to  the  overall  piston  heat  transfer.  Studies 
indicate  that  most  heat  transfer  takes  place  in  the  upper  ring  area 
and  on  the  top  portion  of  oil  cooled  lower  piston  crown  (1-7)*  and, 
therefore,  confirm  that  the  wrist  pin  boss  may  be  treated  as  an  insig¬ 
nificant  factor  in  heat  transfer.  This  reasoning  allows  one  to  treat 
piston  and  valve  heat  transfer  as  a  two  dimensional,  steady  state 
problem. 

If  one  first  examines  the  general  relationships  for  three  dimen¬ 
sional  heat  transfer  in  an  element  as  shown  in  figure  6,  he  can  deter¬ 
mine  the  heat  transfer  (8)  in  each  of  the  three  directions  from 

«x  ■  <#->  <’> 

where  A  =  dy  dx  , 

«y  -  -KA  (-|f)  (2) 

where  A  =  dx  dz  , 
and 

qz  =  -KA  (-Jf)  (3) 

where  A  =  dx  dy  , 

where  T  is  temperature,  K  is  thermal  conductivity, and  q  is  heat  transfer. 
The  three  dimensional  heat  conduction  equation  for  homogeneous  materials 
without  interior  heat  generation  is: 


♦References  are  at  the  end  of  narrative. 


1, 


T  *  T(x,y,z,t) 


Figure  6:  Elemental  Heat  Transfer 
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-«  ,  K  (  ih  t  il  t  A  , 

dx*  dy 

where  t  is  time. 


(4) 


For  steady  state  conditions 


-AI  - 
at 


0  and  equation  4  becomes 


<5> 

dx*  ay  dzd 

The  above  is  known  as  Laplace's  Equation  and  is  often  written 
as  AT  =  0.  If  Laplace's  Equation  is  transformed  to  cylindrical 
coordinates  by  substituting  x  =  r  cos©  ;  y  =  r  sin©  ;  and  z  *  z; 
and  if  cylindrical  synmetry  is  assumed 


the  equation  reduces  to: 


0), 


<)21 

i? 


6  T 

a  r 


(6) 


Since  symmetry  is  justified  in  this  study,  equation  6  becomes 
the  governing  equation  for  pistons  and  exhaust  valves.  Figure  7 
illustrates  this  application  for  a  piston. 

To  solve  the  heat  transfer  problem  for  such  components,  one 
must  solve  Laplace's  Equation  having  certain  boundary  conditions 
known.  A  finite  difference  approach  to  this  solution  has  been  selected 
and  will  follow. 


FINITE  DIFFERENCE  METHOD 

Laplace's  Equation  has  a  unique  solution  for  each  set  of  boundary 
conditions  specified.  However,  analytical  solutions  have  been  worked 
out  only  for  problems  with  certain  special  geometries  and  boundary 
conditions.  The  following  solution  will  be  to  overlay  the  cross  sec¬ 
tion  of  4he  piston  (and  valve)  with  a  grid  of  points  (r,  z)  and  to 
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find  the  temperature  at  each  point  numerically.  A  typical  grid  is 
as  follows  in  figure  8. 

The  numerical  approach  first  approximates  the  partial  deriva¬ 
tives  at  each  point  by  using  a  finite  difference  or 


H 

dr 


(r,z) 


h-j  ,z) 


T 


(r  -  h] 


(7) 


d2T  ,  ,  •  1 

>r2  (r,z>  '  h,2 


[T(r  +  h,. 


z)  +  T(r  -  h1  ,z)  -  2T^r 


,z] 


(8) 


and 


(r,Z)  =  h„2  [T(r’z  +  h2}  +  T(r‘z“h2)  ’  2T(r’z3 


(9) 


where  h^  is  the  horizontal  grid  spacing  and  h 2  is  the  vertical  grid 
spacing. 

Equation  6  now  becomes: 


h  2  jj (r  +  hr  z)  +  T(r  -  h]  ,z)  '  '-T(r,z")j+  rh  [V+hpz)  T(r-h1  ,zj 


4"  [T(r,z  +  h2)  +  T(r,z  -  h2)  '  2T(r,z]  = 


(10) 


For  the  center  line  element  (r  =  0)  a  slightly  different  approxima¬ 
tion  is  needed.  As 


0  ±11 

u’  r  d  r 


d2I 


d  r 


so  the  governing  equation  becomes: 

2  2 

o  .  dT  _ 

2  — J  +  — IT  - 

dr  dz6 


0. 


(H) 
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Theabove  can  be  simplified  if  h-j  =  h ^  and  the  following  substitutions 
are  made: 


T(r  +  h,z)  Ti  +  1 ,  j 


T(r  -  h,  z)=  T,  .  j 


T(r,  z  +  h)  Ti,  j  +  1 


The  following  finite  difference  equations  exist  for  all  points  except 
where  r  =  0: 


0  =  -4  Ti,j  +  Ti  ,j  +  1  +  Ti ,  j  -  1  +  ^  "  7r  Ti-l,j  + 


(  1  +2?}  Ti  +  l,j 


(12) 


For  center  line  points,  (r  =  0),  the  corresponding  finite  difference 
equation  is: 


0  =  -6  T.  .  +  2T.  ,  .  +  2T.  .  .  +  T.  .  .  ,  +  T.  .  ,  ,,,, 

l  ,j  l  -  l,j  i  +  1  ,j  i,  j  +  1  l ,  j  - 1  (13) 


which  by  symmetry  yields: 


(14) 


MODELING  THE  CERAMIC  COATING 

Special  consideration  must  be  given  to  the  points  where  the  ceramic 
coating  meet  the  metallic  piston  (valve)  surface.  Linear  equations  for 
these  pointsare  obtained  by  assuming  a  one  dimensional  flow  across  the 
ceramic  and  using  the  heat  equation: 
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q  =  -KA 


A1 


(15) 


This  is  illustrated  in  figure  9. 

To  obtain  an  equation  involving  T-j ,  T^,  and  T3  one  can  calculate 
the  heat  flowing  through  the  interface  area  (A)  both  from  the 
ceramic  and  the  metal. 

T1  ‘  T2 

q  =  -KCA  (  g— )  (16) 


and 

0  -  -V  (  >  07) 

Equating  these  equations  yields: 


KcA 

-f  (T1  -  V  = 


K  A 

ir  (T2  -  T3^ 


(18) 


or 


T2 


h2  Kc  T1  +  d  T2 

h„  K  +  d  K 
2  c  m 


(19) 


This  gives  a  linear  equation  for  each  grid  point  of  contact  from 
the  ceramic  and  the  metal. 


COMBINED  SOLUTION 

An  equation  now  exists  for  the  temperature  at  each  interior  point 
(r,z)  and  the  surrounding  grid  points.  The  n-equations  involving  n- 
unknowns  have  been  generated.  These  equations  can  be  solved  directly 
by  Gaussian  elimination  or  by  an  iterative  process  like  the  Gauss-Seidel 
method.  In  the  Gauss-Seidel  method  the  temperature  at  one  point  (r,z) 
is  solved  in  terms  of  the  temperatures  at  all  surrounding  grid  points. 

The  result  is  a  temperature  distribution  throughout  the  piston  or 
valve. 


i 


Figure  9:  Ceramic  Metal  Interface 


DETERMINATION  OF  THERMAL  CONDUCTIVITY 

The  thermal  conductivity  of  each  substance  involved  must  be 
accurately  known  over  a  range  of  temperatures  before  any  mathematical 
model  can  be  completed.  The  temperature  dependency  of  aluminum, 
steel,  cast  iron  and  other  metals  is  well  known  and  documented  over  a 
reasonable  range.  Plasma  sprayed  zirconium  oxide,  however,  does  not 
have  a  widely  accepted  conductivity  at  the  temperatures  involved  by 
engine  components. 

Before  completion  of  a  mathematical  model,  therefore,  it  was  nec¬ 
essary  to  determine  such  values.  Samples  of  steel  were  machined  to  a 
standard  sample  size.  Several  of  these  were  plasma  sprayed  with  yttrium 
stabilized  zirconium  oxide  using  the  same  techniques  as  for  engine  parts. 
These  samples,  as  well  as  uncoated  specimens,  were  sent  to  the  Thermo- 
physical  Properties  Research  Laboratory  at  Purdue  University  for  eval¬ 
uation.  The  results  of  this  analysis  were  incorporated  in  the  following 
computer  model  and  are  presented  as  appendix  1. 

THE  COMPUTER  PROGRAM 

The  preceding  section  has  illustrated  a  mathematical  method  to 
describe  the  temperature  at  all  interior  points  for  engine  parts  that 
encounter  two  dimensional  (r,z),  steady  state  heat  transfer. 

A  major  portion  of  this  effort  was  to  incorporate  these  equations 
into  a  computer  program  that  would  give  engineers  a  tool  for  designing 
engines  that  incorporate  ceramic  coated  pistons  and  valves.  The 
following  section  outlines  the  details  of  the  resulting  program. 

Figure  10  details  the  step  by  step  procedure  to  input  component 
details  into  the  program  and  to  acquire  the  resulting  temperature  pro¬ 
files  and  heat  transfer. 

To  implement  the  program,  a  grid  must  be  established  similar  to  the 
one  following  in  figure  11.  The  temperature  at  each  solid  dot  is  needed 
in  the  data. 

The  first  data  card  contains  the  number  of  rows,  the  number  of 
columns,  the  number  of  data  sets,  the  conductivity  of  the  coating,  the 
conductivity  of  the  metal,  and  the  acceptable  iteration  limit.  A  data 


WRITE  OUT 
HEAT  FLOW 


Figure  No.  10:  (continued) 

Program  Flow  Chart 
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set  consists  of  the  boundary  temperatures  for  the  test  of  a  parti¬ 
cular  configuration.  This  information  applies  to  all  the  data  sets 
that  follow.  Each  data  set  is  composed  of  four  sections.  The  first 
section  is  a  single  card  which  contains  the  label  for  the  data  set. 

The  second  section  contains  the  temperatures  along  the  outer  edge. 

For  each  row  there  is  one  card  which  contains  the  edge  temperature, 
and  the  dimensions  to  the  point  from  the  top  and  perimeter  of  the 
component  by  column  number.  The  vertical  dimension  is  recorded  in 
inches  of  depth  while  the  horizontal  dimension  is  stated  only  by 
column  number.  The  format  for  these  cards  is  (F10.3,  F9.3,  12). 

The  temperature  is  recorded  in  degrees  Fahrenheit  (in  agreement  with 
the  units  on  thermal  conductivity)  and  the  distance  is  recorded  in 
inches. 

The  third  section  of  the  data  set  is  composed  of  the  surface 
temperatures  and  column  spacing.  For  each  column  in  the  grid,  a  sur¬ 
face  temperature  and  distance  from  the  center  of  the  piston  valve  is 
given.  The  format  is  (2F10.3).  Again,  the  temperature  is  in  degrees 
Fahrenheit  and  the  radial  location  is  recorded  in  inches. 

The  last  section  of  the  boundary  data  contains  the  temperatures  on 
the  bottom  and  interior  wall  (crankcase  surface)  for  the  piston  and  for 
the  tulip  on  the  valve.  Each  data  card  in  this  section  contains  a  row, 
column,  and  temperature  cf  a  boundary  point  not  covered  in  either  sec¬ 
tion  two  or  three. 

A  data  card  of  this  type  with  row  =  20  and  column  =  20  is  needed 
to  signal  the  end  of  the  data  for  each  set.  The  format  for  these  cards 
is  (2I5.F10.3). 

Note  that  the  distance  between  rows  1  &  2  corresponds  to  the  coating 
thickness.  As  written,  this  program  solves  for  9  rows  and  9  columns. 

This  quantity  can  be  increased  by  changing  the  format  and  programming 
for  the  last  data  card  of  the  fourth  data  section. 

The  output  of  this  program  contains  a  representation  of  the  boundary 
values  inputed  into  the  program,  a  final  temperature  for  each  grid  point, 
and  the  rate  of  heat  flow  vertically  from  the  surface  of  the  piston  (or 
valve)  in  BTU/hr. 


A  card  by  card  listing  of  the  final  version  of  the  program 
follows: 

(see  next  page) 

A  sample  set  of  datacards  used  to  test  the  program  is  listed 
below: 

(see  following  pages) 

These  data  cards  produce  an  output  as  follows: 

(see  following  pages,  table  4-Program  Output). 

The  following  is  a  series  of  two  profiles,  each  of  which  roughly 
outlines  the  geometry  of  the  part  being  tested.  The  top  profile  tab¬ 
ulates  inputed  temperatures  and  the  respective  location  of  each  point. 

The  top  outside  points  has  values  of  r  =  1.63  inches,  z  =  0.00  inches 
and  t  =  280  °F.  Interior  points  are  unknown  at  the  time  of  data  input 
and  are  therefore  printed  as  zero. 

The  calculated  heat  flow  through  the  top  surface  is  written  between 
the  two  profiles  and  in  this  sample  case  equals  366.7  BTU  per  hour. 

The  second  (lower)  profile  displays  the  calculated  temperatures  at 
interior  points  as  well  as  the  input  data  from  the  top  profile.  In 
both  profiles,  the  area  shaded  by  the  letter  "U"  represents  the  crankcase 
environment. 

The  program  gives  one  the  ability  to  calculate  heat  flow  and  the 
temperature  and  distribution  for  any  coated  engine  part  with  two  dimen¬ 
sional  heat  flow. 


OPERATIONAL  TEMPERATURES 


GENERAL 

To  design  engines  incorporating  ceramic  coated  parts,  it  is  desir¬ 
able  not  only  to  have  a  mathematical  model  that  calculates  heat  flow  and 
interior  temperatures,  but  also  one  needs  information  about  actual  sur¬ 
face  temperatures  as  they  relate  to  engine  operating  conditions.  In  addition 
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to  surface  temperatures,  certain  interior  temperatures  are  needed 
to  develop  the  mathematical  model.  For  these  reasons,  a  program 
was  undertaken  to  establish  operational  temperatures  for  a  ceramic 
coated  piston  and  exhaust  valve.  Needed  information  included  temp¬ 
erature  data  at  known  points  for  various  power  levels  for  uncoated 
parts  and  for  parts  with  0.040  inch,  0.070  inch  and  0.100  inch  coat¬ 
ing  of  plasma  sprayed  yttrium  stabilized  zirconium  oxide. 


TEST  FACILITIES 

The  engine  selected  for  acquisition  of  temperature  data  was  the 
same  as  that  used  for  the  previous  performance  data  of  ceramic  coated 
parts.  This  was  a  cooperative  fuel  research  (CFR)  engine  incorporating 
the  low  speed  crankcase  and  the  split  cylinder  head.  The  engine  and 
the  associated  laboratory  are  shown  in  figure  12.  An  engine  schematic 
is  illustrated  in  figure  13. 

Intake  air  was  drawn  through  an  ice  tower  to  control  humidity.  It 
then  passed  through  a  flow  control  system  consisting  of  a  laminar  flow 
element,  a  surge  chamber  and  a  throttle  valve.  By  this  method  a  constant 
humidity,  pressure  and  temperature  of  air  was  available  to  the  engine 
for  all  tests. 

Fuel  for  all  tests  was  commercial  liquified  petroleum  gas  (LPG). 

Lubricating  oil  temperature  was  controllable  by  an  electric  heater 
located  in  the  oil  sump.  Engine  cooling  was  by  the  standard  CFR 
boiling  water  jacket. 

Engine  power  output  was  measured  by  an  indicated  mean  effective 
pressure  (IMEP)  method.  This  technique  incorporates  a  strain  gauge 
pressure  transducer  located  in  the  combustion  chamber,  an  engine  rotation¬ 
al  position  transducer  on  the  crankshaft  and  a  cycle  indicator  located 
on  the  camshaft.  These  three  transducers  provide  the  necessary  data  input 
to  a  PET  2001-8  computer  for  the  solution  of  the  following  equation: 


IMEP 


(20) 


The  average  of  several  of  the  above  solutions  provides  an  average 


ne  Schematic 
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value  of  the  IMEP  for  any  operational  condition.  This  value  is 
recorded  as  engine  power  level. 

Those  familiar  with  temperature  measurement  recognize  the 
problems  in  establishing  accurate  data  at  point  sources  in  ceramic 
substances.  Such  sources  on  the  interior  of  metallic  parts  also  offer 
a  significant  challenge.  Finally,  when  one  locates  these  points  on  a 
reciprocating  piston  (valve)  of  an  operational  engine,  the  full 
complexity  of  the  situation  becomes  apparent. 

The  above  problem  can  be  subdivided  into  two  parts:  1)  what  method 
should  be  used  to  measure  each  point  temperature,  and  2)  how  should 
the  signal  be  transmitted  from  moving  parts  to  the  exterior  of  the  en¬ 
gine? 

Historically,  thermocouples  (including  thin  film  deposited  types), 
deposited  thin  film  resistors,  thermisters  and  Templugs  have  been 
successfully  used  for  general  temperature  measurement. 

Plasma  sprayed  surfaces  are  relatively  porous  which  would  greatly 
reduce  the  probability  of  successful  application  of  thin  film  deposited 
type  devices.  Ceramic  surfaces  that  are  attached  to  a  metal  structure 
make  the  attachment  and  routing  of  leadwires  nearly  impossible. 

In  this  project,  the  ceramic  surface  must  be  reduced  from  0.100  inch 
thickness  to  zero  in  three  stages  by  a  surface  grinding  technique.  Such 
a  process  would  necessarily  destroy  any  surface  device  and  require  the 
reinstallation  and  recalibration  of  a  new  one. 

Templugs  are  a  threaded  metal  insert  that  measure  temperature  by 
hardness  change.  These  devices  measure  only  the  peak  temperature  and 
would  not  be  compatible  with  a  plasma  sprayed  surface. 

Recognizing  the  problems  with  temperature  measurement  on  a  ceramic 
surface  in  the  combustion  chamber  of  a  reciprocating  I-C  engine,  it  was 
decided  that  optical  techniques  offered  the  best  hope  of  success. 

It  is  known  that  all  surfaces  emit  an  infrared  radiation  as  a 
direct  function  of  surface  temperature.  Therefore,  if  this  signal  could 
be  viewed  and  transmitted  outside  the  engine,  a  method  would  have  been 
found  to  instantaneously  measure  surface  temperature  without  touching 
or  affecting  the  surface. 
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The  infrared  signal  transmission  problem  was  solved  by  a  system 
of  water  cooled  probes  that  incorporated  a  glass  fiber  optics  bundle. 

One  bundle  was  designed  so  it  could  view  both  the  top  and  surface  of 
the  piston  and  the  lower  surface  of  the  exhaust  valve.  Both  of  these 
surfaces  were  within  the  chamber.  A  second  bundle  was  used  to  view 
the  exhaust  valve  tulip  (manifold)  surface.  Details  of  these  probes 
and  associated  hardware  follows. 

Figure  14  shows  details  of  an  actual  fiber  optics  bundle  and  the 
water  cooled  probe  into  which  it  is  assembled.  This  probe  is  shown  in 
position  on  the  engine  in  figure  15.  Note  that  the  probe  fits  into  a 
deck  spacer  between  the  engine  cylinder  and  cylinder  head.  This  tech¬ 
nique  allows  the  probe  to  be  moved  in  and  out  to  view  various  radial 
positions  and  to  be  rotated  180°  to  alternately  inspect  the  piston 
and  the  valve  surface. 

Figure  16  displays  details  of  a  second  water  cooled  probe  designed 
to  view  the  exhaust  valve  tulip.  Installation  of  this  probe  into  the 
CFR  cylinder  head  is  seen  in  figure  17. 

The  optical  signal  from  both  of  these  probes  is  transmitted  by 
flexible  fiber  optics  to  a  detector  that  converts  the  I-R  signal  to  an 
electric  voltage.  The  electric  signal  is  read  as  a  temperature  by  an 
analyzer. 

During  the  tests  the  PET  2001-8  computer  seved  as  the  I-R  tempera¬ 
ture  readout  device,  reading  each  probe  temperature  at  a  selected 
engine  position.  For  the  piston,  the  reading  position  was  the  top  dead 
center  (TDC)  position  on  the  exhaust-intake  overlap  event.  The  com¬ 
bustion  chamber  valve  surface  was  read  at  this  same  cyclic  position. 

The  valve  tulip  was  monitored  with  the  valve  closed  at  the  BDC  position 
at  the  end  of  the  intake  stroke. 

SEMIADIABATIC  PISTON 

As  stated  previously,  piston  temperatures  were  required  both  on 
the  combustion  chamber  and  crankcase  surfaces,  as  well  as  on  selected 
points  on  the  interior  of  the  metallic  structure.  Instrumentation 
plans  called  for  the  piston  combustion  chamber  surface  to  be  measured 
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by  an  infrared  probe  while  all  other  temperatures  used  thermocouples. 

The  basic  piston  design  adopted  is  seen  in  cross  section  in 
figure  18  and  in  assembly  in  figure  19.  This  is  a  three  component 
piston  consisting  of  a  steel  crown,  a  phenolic  spacer  and  an  alum¬ 
inum  skirt  in  assembly.  The  crown  is  fastened  to  the  skirt  by  a 
series  of  screws. 

The  piston  crown  is  a  symmetrical  design  (to  insure  two  direction 
heat  flow)  with  access  holes  drilled  for  interior  thermocouple  in¬ 
sertion.  These  thermocouples  were  held  in  place  by  a  ceramic  cement. 
Other  thermocouples  are  attached  to  the  lower  surface  by  peening. 

Figure  20  identifies  the  position  of  the  eleven  thermocouples  used. 
Figure  21  shows  the  piston  crown  with  thermocouples  installed.  Also 
shown,  in  position,  is  the  phenolic  insulator.  Note  that  a  series  of 
holes  have  been  designed  into  each  side  of  the  insulator.  These  holes 
provide  a  screw  type  connection  between  the  fine  wire  (0.015  inch) 
iron  constantan  thermocouples  and  the  stranded  30  gauge  leadwire. 
Thirteen  wires  are  routed  through  the  piston  skirt  shown  in  figure  22 
and  23;  two  of  these  wires  served  as  redundant  commons  while  the  other 
eleven  were  signal  wires  for  each  of  the  eleven  thermocouples. 

THERMOCOUPLE  WIRE  ROUTING 

Leadwire  termination  on  the  piston  was  by  high  temperature  epoxy, 
as  was  the  attachment  to  the  connecting  rod.  Note  in  figure  24  that 
stress  relief  loops  have  been  placed  at  the  junction  between  the  piston 
and  connecting  rod. 

At  the  big  end  of  the  connecting  rod  another  stress  relief  loop 
joined  the  span  between  the  connecting  rod  and  the  swivel  of  a  multi  bar 
support  mechanism.  This  mechanism  was  designed  to  penetrate  the 
crankcase  at  a  sliding-rotating  bearing  and  provide  support  for  the 
thermocouple  leadwires.  An  exterior  view  of  this  device  is  seen  in 
figure  25.  The  base  of  the  multibar  mechanism  is  a  pin  joint  attached 
to  the  engine  baseplate.  The  thermocouple  wires  continue  from  this 
location  through  a  selector  switch  to  computer  and  digital  readout. 

The  exhaust  valve  configuration  without  ceramic  coating  is  shown 


PISTON  CROWN 


PHENOLIC 

“INSULATOR 

SCREW 


PISTON 


Figure 


k 


Figure  25:  Multilink  Support 
Mechanism 
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Figure  24:  Piston-Connecting  Rod 
Assembly 
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in  figure  26.  As  with  the  piston  a  0.100  inch  coating  of  yttrium 
stabilized  zirconium  oxide  was  applied  by  plasma  spray.  Coating 
thickness  was  reduced  to  0.070,  0.040  and  0.000  inches  for  repeated 
tests.  Controlled  valve  seat  contact  was  established  by  hand  grind¬ 
ing,  lapping  and  fitting.  Special  gauges  were  fabricated  to  insure 
repeatabi 1 i ty  of  valve  seat  contact.  The  valve  grid  network  utilized 
in  this  test  is  defined  by  radial  locations  as  viewed  by  the  fiber 
optics  probe  and  as  portrayed  in  figure  27. 

ENGINE  ASSEMBLY 

Figure  28  portrays  the  engine  as  assembled  for  final  test.  Note¬ 
worthy  items  are: 

(A)  Engine  rotational  position  transducers 

(B)  Cycle  position  transducer 

(C)  Multibar  support  mechanism 

(D)  Deck  spacer  with  fiber  optics. 


Engine  test  specifications  follow  in  table  5: 


Table  5 


Bore 

Stroke 

Compression  ratio 

Cooling  water  temperature 

Lubricating  oil  temperature 

Ignition  timing 

RPM 

Fuel 


3.25  inch 
4.50  inch 
6.5  to  1 
210°  F 
120°  F 
10°  BTDC 
630 

Commercial  liquified 
Petroleum  gas 


TEST  PROCEDURE  AND  DATA 

The  data  required  for  ceramic  coating  evaluation  was  as  follows: 

1)  Metal  ceramic  interface  temperature  or  thermocouple 
output  for  numbers  1,  2,  3,  4,  &  9. 

2)  Ring  belt  temperature  or  thermocouple  output  for 
number  10  &  11. 
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3)  Crankcase  surface  temperature  or  thermocouple 
output  for  numbers  5,  6,  7,  &  8. 

4)  Piston  and  valve  combustion  chamber  surface 

or  infrared  output  at  various  radial  positions. 

5)  Valve  tulip  surface  temperature  at  various 
positions . 

The  above  data  values  were  determined  for  ceramic  coating  thick¬ 
ness  of  0.100,  0.070,  0.040,  0.000  inches  at  5  power  levels  correspond¬ 
ing  to  100%,  75%,  50%,  25%  and  motoring.  This  data  is  displayed  in 
the  top  profiles  of  tables  6  &  7,  respectively.  Interpretation  of 
tabulation  has  previously  been  outlined.  Heat  transfer  and  the  cal¬ 
culated  temperature  and  distribution  is  shown  below  the  top  profile. 
Maximum  power  for  this  version  of  the  CFR  engine  corresponded  to  an 
indicated  mean  effective  pressure  of  40.5  psi.  In  each  test,  a  fuel 
flow  was  established  that  produced  this  rated  power.  Fuel  flow  was 
incrementally  reduced  to  produce  the  indicated  steps  of  reduced  power. 
Engine  temperatures  were  allowed  to  stabilize  before  data  was  taken. 

DATA  ANALYSIS 

The  results  of  the  preceding  calculations  are  displayed  in  tables 
8  and  9.  Table  8  displays  heat  transfer  information  based  on  bound¬ 
ary  conditions  using  infrared  data  as  well  as  thermocouples.  Table  9 
is  comparative  data  for  the  piston  only,  based  on  thermocouple  data. 

In  both  cases,  the  data  shows  the  expected  reduction  in  heat  flow 
with  increased  coating  thickness.  The  correlation  between  the  sets 
of  data,  however,  is  not  good.  The  great  difference  in  heat  flow  for 
infrared  and  noninfrared  calculations  leads  one  to  believe  that  the 
accuracy  of  the  surface  temperature  measurement  is  poor. 

ENDURANCE  TEST 

As  previously  stated,  a  part  of  this  project  was  to  establish  the 
endurance  behavior  for  zirconium  oxide  coated  engine  parts.  Actually, 
the  test  described,  herein,  was  the  continuation  of  a  1000  hour  en- 
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Table  7 

Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 


CERAMIC  =  100/1  GOO.  POW6R=  25  per  CENT 


0-00  0.11  0.25  0.38  0.50  3.75  i.QO 


DISTANCE  0.00  0.11  0.25  0.38  0.50 


Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 


1N3D  b3  J  v,  =d3»^cl  ’OOOI/OZ  =  3!Htfc3> 


83 


Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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86 


Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  7  (continued) 
Piston  Temperature  Profiles 


Table  7  (continued) 
Piston  Temperature  Profiles 
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Table  8 

INFRARED  TO  THERMOCOUPLE  HEAT  FLOW 


CERAMIC 

THICKNESS 

INCHES 

PERCENT 

POWER 

HEAT  FLOW 

BTU/HR 

PISTON 

VALVE 

.10 

100 

736 

35 

.07 

100 

888 

61 

.04 

100 

1321 

131 

.00 

100 

16203 

3197 

.10 

75 

706 

27 

.07 

75 

980 

153 

.04 

75 

1619 

233 

o 

o 

75 

17780 

3118 

.10 

50 

703 

78 

.07 

50 

1082 

158 

.04 

50 

1769 

120 

.00 

50 

19763 

3203 

.10 

25 

740 

111 

.07 

25 

1284 

133 

.04 

25 

1780 

220 

.00 

25 

22137 

449 

.10 

MIN 

789 

116 

.07 

MIN 

1399 

160 

.04 

MIN 

1869 

199 

.00 

MIN 

21863 

1579 

i 
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Table  9 

THERMOCOUPLE  TO  THERMOCOUPLE  HEAT  FLOW 


CERAMIC 

THICKNESS 

INCHES 

PERCENT 

POWER 

PISTON 

HEAT  FLOW 

BTU  PER  HOUR 

.10 

100 

2311. 

.07 

100 

2180. 

.04 

100 

2844. 

.00 

100 

2806. 

.10 

75 

1951. 

.07 

75 

2128. 

.04 

75 

2592. 

.00 

75 

2813. 

.10 

50 

2352. 

.07 

50 

1832. 

.04 

50 

2186. 

o 

o 

50 

3646. 

.10 

25 

1632. 

.07 

25 

1502. 

.04 

25 

1936. 

.00 

25 

2696. 

.10 

MIN 

935. 

.07 

MIN 

1159. 

.04 

MIN 

1254. 

.00 

MIN 

1521. 

_ H 


durance  test  performed  under  an  earlier  grant. 

The  purpose  of  this  test  was  to  identify  failure  modes  of 
zirconium  oxide  coatings  in  internal  combustion  engine  combustion 
chambers.  Three  types  of  metal  are  commonly  found  In  engines: 
aluminum,  cast  iron  and  steel.  It  was, therefore,  desirable  to  test 
the  behavior  of  all  three  materials  with  various  coating  thicknesses. 

To  satisfy  the  test  criteria, a  1952  Chevrolet  216  cubic  inch, 

6  cylinder  engine  was  selected.  This  particular  engine  was  chosen 
since  both  cast  iron  and  aluminum  pistons  are  available.  This  engine 
was  remanufactured  and  fitted  with  new  oversize  pistons,  three,  of 
which  were  cast  iron  and  three  were  aluminum.  Prior  to  assembly,  the 
pistons  were  coated  with  a  nominal  0.03,  0.06  and  0.10  inch  coating 
of  yttrium  stabilized  zirconium  oxide.  The  piston  tops  were  machined 
to  reduce  thicknesses  so  that  all  pistons  would  have  the  same  height 
after  coating.  The  valves  and  combustion  chamber  walls  were  also 
coated  with  ceramic  to  the  above  thicknesses.  Figures  29  and  30 
portray  typical  preassembly  component  details. 

The  engine  was  assembled  and  operated  on  a  simulated  driving 
cycle  for  a  first  1000  hours.  This  driving  cycle  consisted  of  23  hours 
plus  55  minutes  at  2000  rpm  at  30%  of  maximum  power,  followed  by  5  min¬ 
utes  at  full  power  also  at  2000  rpm.  No  operational  data,  except  time, 
was  recorded. 

Following  the  first  1000  hours,  the  engine  was  disassembled  and 
inspected.  No  significant  coating  degradation  was  found,  except  on 
the  perimeter  of  the  intake  valves;  here,  most  of  the  ceramic  had 
separated. 

The  second  1000  hours  test  proceeded  after  an  engine  overhaul 
that  included  new  piston  rings  and  bearings. 

The  second  test  was  conducted  with  the  same  simulated  driving 
cycle.  Forty  hours  into  the  second  1000  hour  test,  a  generator  failure 
resulted  in  the  fan  belt  failure.  This,  in  turn,  caused  high  engine 
temperature.  A  coincidental  failure  of  the  safety  switch  allowed  the 
engine  to  run  to  seizure.  Following  this  problem,  the  engine  was 
again  overhauled  and  placed  into  service  with  all  of  the  original  ceram- 


ic  surfaces. 

After  2000  hours  of  exposure,  the  engine  was  disassembled 
and  inspected.  As  before,  no  significant  problems  were  observed 
with  the  coating.  The  piston  surfaces  had  some  heat  cracks  and  some 
feathering  on  the  perimeter;  no  additional  separations  had  occurred, 
however,  since  the  first  1000  hours. 

Figures  31-36  portray  typical  surfaces  after  2000  hours  of  oper¬ 
ation.  Note  that  some  surfaces  seem  abraded  and,  in  some  cases,  this 
condition  goes  to  the  base  metal.  This  abrasion  is  the  result  of 
glass  bead  blasting  that  not  only  removed  carbon  deposits,  but  also, 
the  ceramic  coating.  The  coating  was  secure  and  in  place  prior  to 
blasting. 


SUMMARY  AND  CONCLUSIONS 

The  preceding  document  outlines  details  of  a  study  of  heat 
transfer,  temperature  profiles  and  endurance  for  ceramic  coated 
internal  combustion  engine  parts.  Such  parts  offer  significant  promise 
in  the  advancement  of  engine  technology. 

The  first  objective  of  the  study  was  to  create  a  mathematical 
model  that  could  predict  heat  transfer  and  temperatures  within  coated 
pistons  and  valves.  This  model  was  developed,  and  is  presented,  in 
computer  card  format.  With  this  model,  a  design  engineer  should  be 
able  to  optimize  such  parts  for  any  engine  application. 

The  second  objective  was  to  measure  operational  temperatures  of 
a  piston  and  a  valve  with  various  thicknesses  of  coating  at  various 
power  levels.  To  do  this,  two  rather  novel  techniques  were  used. 
Thermocouples  were  installed  in  moving  crankcase  piston  surfaces  and 
the  signal  was  routed  by  hard  wire  techniques  out  of  the  crankcase. 

This  procedure  was  quite  successful  and  produced  excellent  results. 

The  second  technique  incorporated  fiber  optic  probes  to  measure 
combustion  chamber  surface  temperatures  using  infrared  detectors. 

This  procedure  produced  poor  results  in  spite  of  extensive  calibration 
efforts.  It  is  believed  that  this  method  can  be  made  to  work,  and 
with  additional  development  will  be  an  extremely  valuable  tool  for 


Cylinder  Head  Surfaces 


Figure  35:  Piston  Surfaces 
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engine  research. 

It  is  recommended  that  additional  efforts  be  undertaken  to 
perfect  the  fiber  optics,  infrared  technique. 

The  third  objective  was  the  continuation  of  an  earlier  1000 
hour  endurance  test  of  ceramic  coated  engine  parts.  This  final 
test  added  a  second  1000  hour  of  time  to  the  existing  engine.  This 
effort  proved  without  doubt  that  ceramic  coated  parts  can  endure 
the  torturous  environment  of  I-C  engines  without  significant  degra¬ 
dation. 
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Thermophysical  Properties  of  a  Coating 

INTRODUCTION 

Four  duplicate  samples  of  a  zirconia  coating  on  a  steel  sample  were  furnished 
by  Dr.  Murray  of  Oklahoma  State  University  for  thermal  conductivity  determina¬ 
tions.  The  samples  were  nominally  0,5  inch  diameter  and  0.1  inch  thick  with 
the  steel  comprising  a  little  over  one-half  of  the  thickness.  The  basic  approach 
to  be  used  to  determine  the  thermal  conductivity  "k”  was  to  measure  the  thermal 
diffusivity,  "a",  sjiecific  heat  "C^"  and  density  "d"  since  k  =  aC^d.  The  bulk 
densities  were  determined  from  mass  and  geometry  measurements,  the  specific 
heat  was  measured  using  a  differential  scanning  calorimeter,  and  the  thermal 
diffusivity  "a"  was  measured  using  the  laser  flash  method.  This  technique  for 
obtaining  thermal  conductivity  is  widely  used  and  has  the  advantages  of  simple 
small  sample  size,  and  the  olimation  of  the  heat  flux  and  temperature  gradient 
m.-asurements.  These  later  measurements,  combined  with  extraneous  heat 
flows,  are  the  primary  resons  why  direct  conductivity  measurements  are  diffi¬ 
cult  and  often  yield  poor  results. 

The  Perkin-Klmer  DSC  II  differential  scanning  calorimeter  is  used  for  rou¬ 
tine  determinations  of  specific  heat  up  to  700°  C.  The  differential  heat  required 
to  heat  the  sample  and  a  standard  saphire  specimen  is  measured.  From  the  mass 
of  the  sample  and  the  saphire  and  the  known  specific  heat  cf  saphire,  the  specific 
heat  of  the  unknown  can  be  calculated.  A  continuous  recording  of  microcalories 
versus  temperature  is  obtained  during  the  experiment. 

The  flash  method  in  which  the  front  face  of  a  small  disc-shaped  sample  is 
subjected  to  a  short  laser  burst  and  the  resulting  rear  face  temperature  rise 
is  recorded,  is  used  in  over  80% of  the  present  thermal  diffusivity  measure¬ 
ments  throughout  the  world.  A  highly  developed  apparatus  exists  at  TPRL  and 
we  have  been  involved  in  an  extensive  program  to  evaluate  the  technique  and 
broaden  its  uses.  The  apparatus  consists  of  a  Korad  K2  laser,  a  high  vacuum 
system  including  a  bell  jar  with  windows  for  viewing  the  sample,  a  tantalum 
tube  heater  surrounding  a  sample  holding  assembly,  a  spring-loaded  thermo¬ 
couple,  appropriate  biasing  circuits,  amplifers,  A-D  converters,  crystal 
clocks,  and  a  minicomputer-based  digital  data  acquisition  system  capable  of 
accurately  taking  data  in  the  40  microsecond  and  longer  time  domain.  The 
computer  controls  the  experiment,  collects  (he  data,  calculates  the  results, 
and  compared  the  raw  data  with  the  theoretical  model.  A  typical  visual  display 


is  shown  in  Figure  1.  The  raw  data  are  given  in  Table  1  and  the  computations 
are  given  in  Table  2.  The  vertical  axis  of  Figure  1  is  non-dimensional  tempera¬ 
ture  rise  (the  temperature  rise  at  anytime  divided  by  the  maximum  rise).  The 
horizontal  axis  is  non-dimensional  time  (the  time  from  the  initiation  of  the  laser 
pulse  divided  by  the  time  to  reach  one-half  of  the  maximum  temperature  rise) . 
The  solid  curve  represents  the  theoretical  model  and  the  {Joints  are  the  actual 
data.  Note  the  close  agreement  between  the  experimental  data  and  the 
mathematical  model  ( Figure  l) .  The  calculated  results  given  in  Table  2  were 
obtained  from  the  computer  output  of  Table  1.  The  column  headed  "ALPHA" 
are  the  diffusivity  values  (cm2  sec-1)  calculated  assuming  no  finite  pulse 
effect  at  various  percent  rises  (PER).  The  column  marked  "VALUE"  is  the 
emf  output  corresponding  to  the  percent  rise  and  "'TIME"  represents  the  time 
in  seconds  from  the  initiation  of  the  pulse  until  the  rear  face  temperature  rise 
reached  the  indicated  percent  value.  The  resulting  values  of  a  range  from 
0.00309  to  0.00321  between  10  and  80%  of  the  temperature  rise  (Table  2) 
and  the  best  value  is  0.00320  cm2  sec-1. 

A  special  technique  has  been  developed  to  measure  the  thermal  diffusivity  of 
one  layer  in  intimate  contact  with  a  second  layer.  The  half-time  for  the 
composite  is  measured  and  the  rise  curve  is  compared  to  the  theoretical  model 
in  the  usual  fashion.  In  order  to  use  this  method,  it  is  necessary  to  know  the 
thicknesses,  specific  heats  and  densities  of  the  two  layers  and  the  diffusivity 
of  the  second  layer.  Thus  these  properties  were  measured  for  the  steel  sample 
using  a  blank  furnished  by  Dr.  Murray.  However,  the  rise  curves  for  the 
composite  samples  did  not  follow  the  theoretical  model,  showing  that  inter - 
facial  contact  resistance  was  present.  A  typical. result  is  shown  in  Figure  2. 
Note  the  differences  between  the  theoretical  and  experimental  rise  curves. 

Since  the  coatings  were  reasonably  thick  and  could  be  removed  from  the  sub¬ 
strate,  the  diffusivity  of  the  free-standing  coating  was  measured  directly. 

In  fact  the  visual  display  of  Figure  1  and  the  data  of  Tables  1  and  2  are  from 
the  free-standing  coating  at  room  temperature. 
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RESULTS 

The  thermal  diffusivity  results  for  the  coating  and  for  the  steel  blank  are 
given  in  Table  3.  The  diffusivity  results  for  the  coating  are  plotted  in  Figure 
3.  It  was  noted  that  the  diffusivity  values  were  unstable  with  time  between 
250  and  500°  C  and  consequently  measurements  were  made  during  the  cooling 
cycle.  These  values  were  significantly  less  than  those  attained  during  initial 
heating. 

The  specific  heat  results  for  the  coating  and  for  the  steel  are  given  in  Table 
4  and  plotted  in  Figure  4.  A  large  peak  was  noted  in  the  specific  heat  of  the 
coating  during  the  initial  heating  cycle  (Figure  4).  Thus  the  sample  was  rerun 
and  the  peak  disappeared.  This  behavior  occured  over  the  same  temperature 
range  in  which  the  diffusivity  values  were  unstable  (Figure  3)  and  the  two 
phenomena  are  undoubtedly  related. 

Thermal  conductivity  values  are  calculated  in  Table  5  and  plotted  in  Figure 
5.  No  calculations  were  made  between  300  and  500°  C  for  the  initial  heating 
data  as  the  specific  heat  and  diffusivity  values  are  certainly  time -temperature 
dependent.  However  the  cooling  curve  data  (diffusivity)  and  second  run 
data  (specific  heat)  are  representative  of  this  material  after  temperature 
cycling  to  500°  C  ( 1000°  F) . 

With  the  knowledge  of  the  thermal  properties  of  both  layers,  it  is  possible 
to  calculate  the  contact  conductance  for  the  two-layer  case.  These  values 
ranged  from  2.4  to  4.2  W  cm-2  K-1  at  room  temperature  for  the  three  samples 
which  were  not  sectioned. 

The  conductance  values  at  higher  temperatures  or  the  changes  caused  by 
temperature  cycling  were  not  determined,  although  it  would  be  possible  to 
do  so. 

DISCUSSION 

The  changes  in  specific  heat  and  thermal  diffusivity  (and  consequently 
conductivity)  noted  during  initial  heating  are  undoubtedly  due  to  the  removal  of 
a  contaminant.  We  have  often  observed  similar  phenomena  on  other  samples 
of  materials  which  were  significantly  less  than  theoretical  density.  The  sensi¬ 
tivity  and  accuracy  of  advanced  computer-aided  instrumentation  and  methods 
makes  such  changes  readily  discernable.  However,  unless  one  is  interested  in 
initial  heating  conditions,  one  should  use  the  values  generated  during  the  cooling 
or  second  run  curves. 


TABLE  1 


Experimental  Data 


OKLA. 

RUN  #  1  ON  7/29/71. 

DURATION:  2.  200000  SEC.  ; 

SHORT  BASELINE; 


STATE  UNIV.  COATING  AIR 
SAMPLE  THICKNESS:  O 
0. 10000  SEC.  BETWEEN 
NO' WILDNESS  &  SMOOTHING 


VOLTAGE 
O.  47363159 
O.  50043699 
0.  51269400 
O.  55903060 
O.  86181420 
1.  27196937 

1.  716304  19 

2.  20214  277 
2.  63915337 

2.  95653536 
3. 29509001 

3  60833919 

3.  81102539 
4. 0087 7873 

4.  19676661 
4.  37254734 
4. 44823078 

4  62157013 
4.  65086698 
4.  71922613 
4.  83335474 

4  89500701 

4.  94627632 

5.  03172539 
5.  03172539 
5.  02928398 
5.  09520130 

5  14158837 
5.  12449360 
5.  17083517 
5.  19885733 
5.  13553361 
5.  21838900 
5  84163500 
5.  19285  733 
5.  21971322 
5.  29733800 
5.  26.121698 
5.  85145135 
5.  32957620 


TIME 

0. 

44300. 
99610. 
154920. 
210230. 
265530. 
320830. 
376150. 
431460. 
486760. 
542060. 
577370. 
652630. 
700000. 
763300. 
81.3600. 
873910. 
929220. 
784530. 
1039830. 
1095130. 
1150440. 
1205750. 
1261050. 
1316360. 
1371660. 
1426970. 
14  82230. 
1537590. 
1592890. 
1648200. 
1703510. 
1758020. 
1014120. 
10694  20. 
1924730. 
1930050. 
2035360. 
2090660. 
2145960. 


INITIAL  VOLTS 
0.  43827999 
0.  48339720 
0.  512694  00 
O.  55175639 
0.  48827999 
O.  507811  19 
0.  47363159 
0.  43827999 
O.  54199080 
0.  522.4  5959 
O.  55175639 
0.  49316280 
O.  51269400 
0.  51269400 
0.  59326019 
0.  58105319 
0.  54179080 
O.  51269400 
0.  52245959 
0.  50292839 
0.  49316280 
0.  55175639 
0.  54199090 
0.  43339720 
0.  49316280 
O.  51757679 
O.  53222519 
O.  50292839 
O.  53222519 
0.  52734240 
0.  49316280 
0.  48827999 
0.  50292839 
0.  50701  119 


23C 

1046  CM; 

TRAILING  READINGS 
PERFORMED; 


FINAL  VOLTS 
5.  31004503 
5.  34666597 
5.  33690042 
5.  37596277 
5.  44920481 
5.  40525957^ 
5.  4  1502520 
5.  39549395 
5.  4  5408762 
5.  44920481 
5.  4  7361880 
5.  46385318 
5.  47361880 
5.  48338435 
5.  43455637 
5.  46873599 
5.  47361880 
5.  50779834 
5.  47361880 
5.  44432200 
5.  45408762 
5.  47361880  ’ 
5.  4  34  55637 
5.  47361880 
5.  43455637 
5.  43455637 
5.  41502520 
5.  38572S32 
5.  41502520 
5.  4  5408762 
5.  43455637 
5.  37596277 
5.  39549395 
5.  42967364 


TABLE  2 


Calculation  of  Thermal  Diffusivity  Values 
from  Experimental  Data 


OKLA.  STATE  UN1V.  COATING 


VALUE 

EASEL  I NE:  O.  51312470 

HALFMAXI:  2.92566239 

MAXIMUM:  5.31561374 


ALPHA 

PER 

0.  00315968 

10.  0 

0. 00320300 

20.  0 

0. 00321028 

25.  0 

0.  00320497 

30.  0 

0.  00321118 

33.  3 

0.  00320276 

40.  0 

0.  00319483 

50.  0 

0. 00318505 

60.  0 

0  00317127 

66.  7 

0.  00310430 

70.  0 

0.  00309031 

75.  0 

0.  00312028 

SO.  0 

0  00300031 

90  O 

FLASH.  FT 

DATE:  7/29/71 

\ 

AIR  23C 

TIME (SEC) 

0. 40661017 

0.  47707939 

0.  40919011 

VALUE 

TIME (SEC) 

0.  99337363 

0.  22891539 

1.  47362291 

0.  28780382 

1. 71374678 

0.  31602185 

1.  95387184 

0.  34  552216 

2.  11395442 

0.  36448961 

2.  434  12017 

0.  40638661 

2.  91436910 

0.  47528982 

3.  39461803 

0.  55730617 

3.  71478438 

0.  62469685 

3.  87486696 

0.  67626202 

4. 11499202 

0.  74320042 

4. 3551 1648 

0.  81770313 

4  33536422 

1  10555171 
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TABLE  3 

Thermal  Diffusivity  Results 


Sample 

(No.) 

Temp. 

(°C) 

Temp 

(°F) 

Steel 

23. 

73 

Steel 

93 

199 

Steel 

204 

399 

Steel 

315 

599 

Steel 

426 

799 

Steel 

537 

999 

Coating 

23 

73 

Coating 

93 

199 

Coating 

204 

399 

Coating 

315 

599 

Coating 

426 

799 

Coating 

537 

999 

Coating 

426 

799 

Coating 

315 

599 

Coating 

205 

399 

Coating 

93 

199 

Coating 

23 

73 

Diffusivity 
(cm2  sec"4) 

0. 183 
0. 163 
0.137 
0.113 
0.  094 
0.078 


0. 00320 
0.00305 
0. 00290 
0. 00274 
0. 00257 
0. 00241 
0. 00244 
0. 00251 
0. 00258 
0. 00273 
0. 00288 
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TABLE  4 


Specific  Heat  Results 


Temp. 

Temp. 

C  oating 

Steel 

(°K) 

<°C) 

Run  1 

Run  2 

350 

77 

0.51905 

0.  52145 

0.47724 

375 

102 

0. 53324 

0. 54027 

0.48175 

400 

127 

0.54356 

0.54995 

0.49227 

425 

152 

0.55271 

0. 56119 

0.50173 

450 

177 

0.56190 

0.57115 

0. 50846 

475 

202 

0.56759 

0.57716 

0.51875 

500 

227 

0.57740 

0. 58447 

0.52755 

525 

252 

0.58942 

0. 59532 

0.53683 

550 

277 

0.60218 

0. 59873 

0. 54504 

575 

302 

0,61699 

0. 60142 

0. 55349 

600 

327 

0. 63429 

0. 60287 

0. 56312 

625 

352 

0. 65684 

0.60469 

0.57390 

650 

377 

0.68326 

0.60510 

0. 58512 

675 

402 

0.  71580 

0. 60443 

0. 59739 

700 

427 

0. 75878 

0. 60494 

0.60973 

725 

452 

0.81170 

0. 60692 

0.  62316 

750 

477 

0. 85459 

0. 60988 

0.63734 

775 

502 

0.86087 

0.61402 

0.65170 

800 

527 

0. 86361 

0. 62104 

0. 66687 

825 

552 

0.  85323 

0. 63070 

0. 68309 

850 

577 

0.80381 

0.64010 

0. 70189 

875 

602 

0. 75636 

0.64876 

0.  72623 

900 

627 

0.72906 

0.65818 

0. 75374 

925 

652 

0.  71783 

0.  66983 

0.78855 

<;  s 
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Figure  1.  Comparison  of  Experimental 
Rise  Curve  with  Theoretical  Model 
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